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Abstract 

We study the exclusive decay of J/^f, T and B c into a lepton pair combined with two 
pions in the two kinematic regions. One is specified by the two pions having large momenta, 
but a small invariant mass. The other is specified by the two pions having small momenta. 
In both cases we find that in the heavy quark limit the decay amplitude takes a factorized 
form, in which the nonperturbative effect related to heavy meson is represented by a NRQCD 
matrix element. The nonperturbative effects related to the two pions are represented by some 
universal functions characterizing the conversion of gluons into the pions. Using models for 
these universal functions and chiral perturbative theory we are able to obtain numerical 
predictions for the decay widths. Our numerical results show that the decay of J/ip is at 
order of 10 -5 with reasonable cuts and can be observed at BES II and the proposed BES III 
and CLEO-C. For other decays the branching ratio may be too small to be measured. 
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1 Introduction 



50 x 10 6 J/ip events have been collected with the upgraded Bejing Spectrometer (BES II) at Beijing 
Electron Positron Collider (BEPC), and several billions J/ip events will be collected with the 
proposed BES III at BEPC II and CLEO-C at modified Cornell Electron Storage Ring (CESR)H, 

Furthermore, about 4 fb _1 bb resonance data are planned to be taken at CLEO III in the year 
prior to conversion to low energy operation (CLEO-C) [0]. With these data samples various decay 
modes of J/tp and bb resonances can be studied with high statistics. In this paper we propose 
to study the exclusive decay of 1 quarkonia and B c into a lepton pair and a pion pair. We 
consider two limited cases in the kinematic region. One is specified by the pion pair having a 
large total momentum and a small invariant mass. In this case, the pions are hard. The other is 
specified that the pion pair having a small momentum, i.e., the pions are soft. In these decays 
the two pion system is produced through conversion of gluons into the two pions. Because of 
isospin symmetry conversion of gluons into one pion is highly suppressed. In a two-pion system 
the two pions can be in an isospin singlet, the conversion is allowed. Hence these decays will 
provide valuable information how unobservable gluons, as dynamical degrees of freedoms of QCD, 
are converted into observable hadrons. 

In the case that two-pion system has an invariant mass m n7T which is much smaller than the 
mass of heavy meson and has a large total momentum, the decay amplitude takes a factorized 
form in the heavy quark limit, in which the nonperturbative effect related to heavy meson is 
represented by a non-relativistic QCD (NRQCD) matrix element 0, and that related to the two 
pions is represented by a distribution amplitude of two gluons in the isoscalar pion pair which is 
defined with twist- 2 operators. The two gluons are hard in the kinematical region, their emission 
rate can be calculated with perturbative QCD. The same distribution amplitude also appears in 
the predictions for productions of two pions in exclusive processes 7 + 7* — > 7r+7rj|, |5|, f|, j*+h — > 
h + 71 + 71 [0, H , and the radiative decay of 1 heavy quarkonium ||, where the amplitudes can be 
factorized in certain kinematic region. Besides these processes, the decays studied here will provide 
another way to study the nonperturbative mechanism of how gluons, which are fundamental 
dynamical freedoms of QCD, are transmitted into the two pions. Furthermore, for 7 + 7* — > 7T+/T, 
at the tree-level, only the distribution amplitude of quark appears in the scattering amplitude, 
the distribution amplitude of gluon appears at loop levels or through evolution of distribution 
amplitudes]!, B 0' while for 7* + h — > h + n + 71, at the tree-level, the distribution amplitude of 
quark as well as that of gluon contribute to the scattering amplitude, but the produced charged 
pion pair is dominantly in an isospin 1=1 state |7], ||. This may make it difficult to extract 
the distribution amplitude of gluon from corresponding experimental data, because the two pions 
produced through gluon conversion are in a I = state. For the decays studied in this paper and 
the radiative decay of 1 heavy quarknonium to two pions|§, only the distribution amplitude of 
gluon appears at the tree-level and the produced two pions are dominantly in a I = state. This 
makes the extraction of gluon content of a two-pion system relatively easier in experiment. Of 
course, comparing with the radiative decay of 1 heavy quarkonium into two pions in the same 
kinematic region, the leptonic decay of 1 heavy quarkonium to two pions is suppressed by the 
fine structure constant a, but the final state in the latter case is more clearer and can be detected 
with higher efficiency. With the model for the distribution amplitude of gluon, developed in || |] . 
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we obtain numerical predictions for the branch ratio of the decay in the considered kinematic 
region. Our results show that the decay mode of J/ip in the considered kinematic region is 
certainly observable at BES II and the proposed BES III and CLEO-C. For other decays the 
branching ratios may be too small to be measured. 

In the case with tow soft pions, it has been shown that the decay amplitude of J/^f and of T also 



takes a factorized form in the heavy quark limit [ 10 1. In the decay amplitude, the nonperturbative 



effect related to heavy quarkonium and that related to pion pair can be separated, the former is 
still represented by a NRQCD matrix element, while the later is represented by a matrix element 
of a correlator of electric chromofields which characterizes soft gluons transition into the pion 
pair. This result is nonperturbative. For B c decay one can generalize the approach and obtain 
a factorized form for the decay amplitude, where the same correlator appears. Since the matrix 
element of the correlator of electric chromofields between the vacuum state and the two-pion state 



is unknown, no numerical prediction of the decay is given in 10 1. In this paper, we develop a model 



for the matrix element of the correlator of electric chromofields and give numerical predictions for 
the leptonic decays J/ip,T(lS) and B c into a soft pion pair. Numerical results are obtained in 
the considered region of kinematics and show that the decay mode of J/ip is observable at BES 
II and at the proposed BES III and CLEO-C. 

The paper is organized as the following: In Sec. 2 we study the decays of J/ip and T(IS') into 
two hard pions combined with a lepton pair, where we give a detailed derivation of the factorized 
amplitude of the decay. Numerical results for the decays are presented. In Sec. 3 the decays 
of J/ip and T(15) into two soft pions combined with a lepton pair are studied, a model for the 
matrix element of the correlator of electric chromofields is developed, and numerical results for 
the decay are given. Sec. 4 is devoted to the study of the decays of B c . We summarize our work 
in Sec. 5. 

In this paper, we take nonrelativistic normalization for the heavy meson states and for heavy 
quarks, and we take the pion pair to be of a tt + and a 7r~. Using isospin symmetry one can easily 
obtain results for a pair of 7r°'s. 



2 Leptonic decays of J/ip and T(15) to two hard pions 

We study the exclusive decay in the rest frame of J/ip: 

j/iP(p) - i + ( Pl ) + r(pa) + 7i + {K + ) + tt-(K-), (1) 

where I = e, /i, the momenta are indicated in the brackets. We denote k = k n + + k n ~ , q = p\ + p 2 
and = k 2 . We consider the kinematic region where |k| 3> m^n an d <C M^. At leading 
order of QED, the S-matrix element for the decay is 

(f\S\i) = -iQ c e 2 L^ • 1 J d 4 ze^ z (n + n-\c(z)rc(z)\J/iP), (2) 

where Q c is the electric charge of c-quark in unit of e, c(z) is the Dirac field for c-quark, and 

L ll = u(p 2 h fl v(p 1 ). (3) 
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At leading order of QCD, two gluons are emitted by the c- or c-quark, and these two gluons will 
be transmitted into the two pions. Using Wick theorem we obtain: 

(f\S\i) = l\s ab Q^9 2 s L,-j 2 

x j d 4 z d 4 y (fx fx x <fyi e i{q - z+k2 ' y) 

X (Olc^C^lJ/^) (7T + 7T-|G^(xX(0)|0) 

x [5\x - x x )8\z - yi)^S F (x - y)^S F {y - z)Y + • • •]*, (4) 

where is the momentum of one of emitted gluons, Sf{x — y) is the Feynman propagator of 
c-quark, the dots in the square bracket denotes other five terms. In the limit of m c — ► oo, a c- 
or c-quark moves with a small velocity v, this fact enables us to describe nonperturbative effect 
related to J/ip by NRQCD0. For the matrix element (0\cj(xi)Ci(yi)\J/if)) the expansion in v can 
be performed with the result: 

(OlcfCsiMyOIJM = -1 (P+ V P-h (0| X VV| J If) e~^ + y^ + 0(v 2 ), (5) 

where xK^P) * s the NRQCD field for c(c) quark, cr(£ = 1, 2, 3) is the Pauli matrix, and 

ft= = |(1±7°) 
p = (m c , 0,0,0). (6) 

The matrix (O|xV , 0|"V'0) is proportional to the polarization vector e e (J/ip) at the considered 
order. In this paper, we neglect the contribution from higher orders in v, the momentum of J/ip 
is then approximated by 2p. It should be noted that effects at higher order of v can be added 
with the expansion in (H). 




Figure 1: One of the Feynman diagrams for the exclusive decay of J/ip into lepton pair and two 
pions. 

Using (RI) we can write the S-matrix element as 

(f\S\i) = ^Q c e 2 g 2 s (2vr) 4 5\2 P - k - q)L, • I (0| X V^| J/f 

x / r^(Mi), (7) 
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and 

T^(k, ki) = J d 4 x e~ iklX (7r + 7r-|G a '' i (x)G a,,/ (0)|0>, (8) 

where H £fltJ-lUl (p,k,ki) is the amplitude for a cc pair emitting a virtual photon and two gluons, 
and this can be calculated with perturbative QCD. The contributions in (|7|) can be represented 
by Feynman diagrams. One of them is given in Fig.|I|, where the kinematic variables are also 
indicated. The nonperturbative object V iu {k, k\) describes how two gluons are converted into the 
two pions. 

If the two pion system have a large momentum and a small invariant mass, a twist expansion 
for the nonperturbative object r M!/ (fc, k{) can be performed. For convenience we will work in the 
light-cone coordinate system, in which the components of k are given by 

fc M = (k + , k~, 0), k + = (k° + k 3 )/V2, k~ = (k° - k 3 )/V2. (9) 
In the light-cone coordinate system we introduce two light cone vectors and a tensor: 

= (0,1,0,0), ^ = (1,0,0,0), 

= - n »n u -n"*,", (10) 

and we take the gauge 

n-G(x) = 0. (11) 




0.3 0.4 0.5 0.6 0.7 



m nm ( GeV ) 

Figure 2: The differential decay branching ratio of J/ip, d BR(J/^, H)/dm n7T as a function of m w 
in unit of 10 _6 GeV _1 with the cuts given in the text. 

The x-dependence of the matrix element in (|8|) is controlled by different scales. The in- 
dependence is controlled by k + , which is large in the kinematic region we considered, while the 
x + — and XT-dependence are controlled by the scale k~ and Aqcd, which are small in comparison 
with k + . With this observation we can expand the matrix element in x + and in xt- The resulted 
twist expansion of the Fourier transferred matrix element T^ l/ {k, k\) is a collinear expansion in k~ 
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and kx ~ ^qcd- Hence the expansion parameters of Y lLV {k, hi) are k~/k + and Aqcb//c + , with 
k~/k + < 0.10 and Aqcb/A; + « 0.20 for J/-0 in the kinematic region considered. At the leading 
order only twist-2 operators contributes to the matrix element. We will neglect higher orders in 
the expansion, i.e., we only keep contributions of twist-2 operators. Then we obtain: 



r""(fc,fci) = (2vr) 4 5(k7) 5 2 {k 



IT 



k + X\(l — x\) 



d T v $ 



G, 



x 1 ,(,m 7[ 



(12) 



with 



$ G (xi,C,m, 



1 

k+ 
k\ ■ n 

k ■ n 



dx 

c 



x (7i + n-\G a ' + ^(x~n)G a ' + {0)\0), 



n 



n 



(13) 



$ (xi,£,m W7T ) is the gluonic distribution amplitude which describe how a pion pair with helicity 
A = is produced by two collinear gluons; this represents a nonperturbative effect and can only be 
calculated with nonperturbative methods or extracted from experiment. As it stands, it is gauge 
invariant in the gauge n ■ G(x) =0. In other gauges we need to supply a Wilson line operator to 
make it gauge invariant. With (|12|) it is straightforward to obtain the S'-matrix element at the 
tree-level in our approximation: 



(f\S\i) 



— i 

24 



1 



Q c e^(2vr) 4 6\2p - k - q)IS ■ ^(0\ X ^\ J/V) 



o 



1 $ G (xi,C,m 7r7r ) 
ax\- 



x x (l - x Xj 



\d^ 



(p, k, kj 



with 



-d^ 1Vl ■ H, 



(p, k, ki 



16 

Ml 



16 



(14) 



(15) 



where is the mass of J/tp. In ( |15|) we have neglected the mass m W7r , since the effect of m n7T 
should be combined with effects of twist-4 operators as a correction to the above result. In 
( |14] ) the nonperturbative effect related to J ftp and that to the two-pion system are separated, 
the former is represented by a NRQCD matrix element, while the later is represented by the 
distribution amplitude of two gluons in the isoscalar pion pair, which is defined in ( p^ ) in the 
gauge n ■ G(x) = 0. 

The kinematics of the decay can be fully described by five variables as in K e ^ decay [ffT|: 

1. m^ n , the invariant mass squared of the pion pair; 

2. q 2 = (pi +P2) 2 , the invariant mass squared of the lepton pair; 

3. 6^, the polar angle of the tt + in the rest frame of the pion pair with respect to the moving 
direction of the pion pair in the J/ip rest frame; 

4. 61, the polar angle of the Z + in the rest frame of lepton pair with respect of the moving 
direction of the lepton pair in the J/ip rest frame; 
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q ( GeV ) 



Figure 3: The differential decay branching ratio of J/ip, d BR(J/ip, H)/dq 2 as a function of q 2 in 
unit of 10 _6 GeV -2 with the cuts given in text. 



5. <p, the azimuthal angle between the two planes in which the pion pair and the lepton pair 
lies respectively. 



In terms of these variables, the differential decay width can be written as 

1 7T 2 Ijfcl „ „ = 



dT 



(2tt) 8 32 Mjp 
where (3 and Pi are defined as: 



P Pi y^}M\ 2 dq 2 dm 2 n d cos 9 n d cos 9\ 



(16) 



P 



\ 



1 



4 ml 



mi 



Pi 



\ 



4 mf 



(17) 



X]|M| 2 is the absolute squared matrix element of the decay, summed over final state spins and 
averaged over initial state spin. From ([14]) and (|15|) , we have 

512g 2 [(M 2 + g 2 ) + (M 2 -g 2 )cos 2 ^] 



^Q c v^-i-i(oixWiw)r 



x ■ 



x 



3M 2 (M 2 -q 

1 $ G (a;i,C,m 7r ^) 

dxx — 

o xAl — Xi) 



2\2 



which is independent of the azimuthal angle 0, the spin average for J / ip is implied in the squared 
matrix element. 

To give numerical predictions, the nonperturbative inputs, the NRQCD matrix element and 
the distribution amplitude of two gluons in the isoscalar pion pair, are needed. The NRQCD 
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Figure 4: The decay branching ratio of J/ip, BK(J/ip, H, M 2 ) as a function of M 2 with M 2 > q 2 
in unit of 10 -6 . The other cuts are the same. 

matrix element is related to the wave-function of J/ip in potential models and can be estimated 
with these models. It can also be calculated with lattice QCD or extracted from experiment. In 
this paper, we use leptonic decay of J/ip to determine the NRQCD matrix element, i.e., 




m ( GeV ) 

Figure 5: The differential decay branching ratio of T(IS'), d BR(T, H)/dm n7T as a function of 
in unit of 10 _7 GeV _1 with the cuts given in the text. 

The distribution amplitude Q G (xi, (, m nn ) is not determined with experiment, a detailed study 
of $ G (a;i, (, m n7T ) can be found in ^ |6], §]. For our numerical prediction we will use their re- 
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suits for Q G (xi, C, rn-xTr), in which asymptotic form of $ G (xi, £, m W7r ) is taken as an Ansatz for 
$ (xi,(, m mr)- It should be noted that the renormalization scale \i should be taken as in 
our case. Because it is not very large, the actual shape of § G (x 1 , £, m nn ) may look dramatically 
different than that of the asymptotic form. Keeping this in mind we take the form <& G (xi, (, m nn ) 
as that given in H: 



$ G (a;i, C, m™) = -60 M G x\ (1 - x x ) 2 
where ( is related to 9 n and m nn by 



3C — 3 2 (3 2 

— ^ — fo(m n7r ) Po(cos^) — — f 2 {m 7Tn ) P 2 {cos6 n ) 



(3 cos6 n = 2 C - 1, 

C is a constant and takes C = 1 + b m 2 + O(m^) with b ~ — 1.7GeV~ 2 
by gluon fragmentation into a single pion, its asymptotic value is 



M. 



G 



4 Ct 



N f + 4 C F 



(20) 



(21) 

M G is determined 



(22) 



f 0(171^) and f 2(171^) are the Omnes functions for I = s- and d-wave hit scattering, respectively. 
The Omnes function f2(m n7T ) is dominated by the / 2 (1270) resonance resulting a peak at = 
1.275GeV, while the Omnes function fo(m n7T ) in the relevant m n7T region we studied (m nn < 



0.70GeV) can be calculated by the chiral perturbative theory, the result is \\L2 



1 + 



mt 



+ 



2ml 



mt 



1927T 2 /. 2 32vr 2 / 2 
where /„- = 93 MeV is the pion decay constant. 



/31n 



1 + 6 



+ 2 + m(3 



(23) 




q ( GeV ) 



Figure 6: The differential decay branching ratio of Y(IS'), d BR(T, H)/dq 2 as a function of q 2 in 
unit of 10~ 7 GeV~ 2 with the cuts given in the text. 
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Cuts must be used to select the kinematic region where the two pion system has a large 
momentum and a small invariant mass. We use the cuts: k + > 10 k~, k° + |k| > 2.0GeV and 
2m w < m ww < 0.70GeV, this corresponds to q 2 < 2.5GeV 2 for J/ip and q 2 < Q7GeV 2 for T(15). 

With these results we are able to predict the decay branching ratio in the considered re- 
gion. The quark masses are take as m c = and mj = \M^. a s {2m c ) = 0.31 for J/ip, 
a s (2mb) = 0.21 for T(1S). All other parameters needed are taken from ]14|]. The m nn distribution 
of J/ip decay, integrated over | cos^| < 1.0, | cos^| < 1.0, < < 2n, and Amf < q 2 < 2.5GeV 2 , 
denoted as d BR(J/ip, H)/dm n7T , is shown in Fig. |2|. The q 2 distribution of J/ip decay, inte- 
grated over |cos07r| < 1.0, |cos#z| < 1.0,0 < <f> < 2n, and 2m.,, < < 0.70GeV, denoted as 
d BR(J/ilj,H)/dq 2 , is shown in Fig. |^. The q 2 distribution decreases rapidly as q 2 increases, this 
is mainly due to the q~ 2 factor of the photon propagator. This behavior is shown in another way 
in Fig. [|, where the decay branching ratio as a function of the cut M 2 with q 2 < M 2 . We see from 
this figure, for M 2 = 10" 2 , 10" 1 , 10°GeV 2 , the corresponding contributions are 59%, 75%, 92% to 
the branching ratio in the whole region considered, respectively. Integrating over the kinematic 
region we considered, the decay branching ratio is 1.0 x 10~ 5 for J/ip, among which the s-wave and 
d-wave contributions are 9.2 x 10 -6 and 8.6 x 10~ 7 respectively, i.e., the d-wave contribution is 
suppressed in the kinematic region here. The results indicate that this decay mode as well as the 
distributions of m nn and q 2 can be observed at BES II, and at the proposed BES III and CLEO-C. 
In above numerical calculations, the renormalization scale of the effective QCD coupling is taken 
to be J/ip mass with Aqq D = 280MeV[14], i.e., a s (2m c ) = 0.31. If this scale is taken to be m c , 
the corresponding decay branching ratio of J / ip in the considered kinematic region increases by a 
factor of 2.1 by using (18). Hence, the decay branching ratio is conservatively predicted. 




1 i — — — — — i 

10 2 10 1 1O 10 1 10 2 

M 2 ( GeV 2 ) 

Figure 7: The decay branching ratio of T(IS'), BR(T, H, M 2 ) as a function of M 2 with M 2 > q 2 
in unit of 10~ 7 . The other cuts are the same. 

The corresponding differential decay branching ratios for T(IS') decay are shown in Figs. |5|- 
[7[ The kinematic region we studied for T(IS') decay is | cos 6*^ | < 1.0, | cos ^ | < 1.0,0 < < 
2n,2m n < m wn < 0.70GeV, and Am 2 < q 2 < 67GeV 2 . It is interesting to observe that in Fig. 
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U there is a turn-over near q 2 = 40GeV 2 . Similar to J/tp case, the dominant contribution to 
the decay of T(1S) comes also from small q 2 region. For M 2 = 10 -2 , 10 _1 , 10°, and lO^^GeV 2 , 
the corresponding contributions are 44%, 56%, 68%, and 81% to the branching ratio in the whole 
region considered, respectively. Integrating over the kinematic region we considered, the decay 
branching ratio is 3.0 x 10~ 7 , among which the s-wave and d-wave contributions are 2.8 x 10~ 7 
and 2.6 x 10 -8 respectively, i.e., the d-wave contribution is also suppressed in the kinematic region 
here. For the decay of T one may allow to be larger than that in the decay of J/ip, i.e., 
0.7GeV, because the phase space is large. With a large upper cut for m^, the branching ratio can 
become large. But fo(m n7T ) determined with chiral perturbation theory may become unreliable 
for large m n7T . 

It should be noted that the two-pion state is produced with the helicity A = at the level of 
leading twist. It can be a mixture of states with different angular momenta L. Because of parity 
conservation, isospin symmetry and Bose-Einstein statistics, L can only be even. Our numerical 
results show that the state is mainly in a s-wave state. At levels of higher twist it is possible that 
the two-pion state is produced with A ^ 0. Following the analysis for the radiative decay of T 
into /2 (1270) one can show that the two-pion state can be produced with |A| = 1 and 2 at 
order of twist 3 and of twist 4, respectively. 



3 Leptonic decays of J ftp and T(lS) combined with two 
soft pions 

In this section, we study the leptonic decays of J/ip and T(1S) combined with two soft pions. 
In contrast to the decays studied in the last section the gluons, which are emitted by the heavy 
quarks in the quarkonium and are converted into the pions, are soft. The emission of soft gluons 
can be studied by employing an expansion in the inverse of the heavy quark mass itlq. It is 



shown in |T0| that at leading order of the expansion the decay amplitude in this kinematic region 
can be factorized into three parts: the first part is a NRQCD matrix element representing the 
bound-state effect of heavy quarkonium, the second part is a matrix element of a correlator of 
electric chromofields, which indicates the nonperturbative effect of the soft gluons converted into 
the soft pion pair, the third part consists of some coefficients. It should be emphasized that the 
results can be derived without using perturbative QCD. In this section we present a model for 
the matrix element of the correlator of electric chromofields, and give numerical predictions. The 
S-matrix for the J/ip decay is[ffD] 



(f\S\i) = i"-Q c e l W 5\2p - k - q)L, ■ ^-(0| X V^| J/il>) 
3 q l 



r 1 l 

T n Jk) = dr : rh(r,k), (24) 

nwK ' J l + r-i0+ 1-t-z0+ V ' K J 

where the momenta are denoted as the same in the last section. For soft pions we have |k| <C rrig 
and <C nig. h(r, k) is the distribution amplitude for the soft gluons converted into two soft 
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pions. It is defined as 



h{r, k) 



9l 
2tt 



dte lTk0t (-K + 7r-\E a {t,0) 



P exp 



dx°G°' c (x , V 



ab 



B\-t,0)\0), 



(25) 

where P denotes path-ordering and r c is the generator of SU(3) in adjoint representation, {T c ) a b = 
—ifabc- Because of the energy conservation h(r,k) = if |r| > 1. The term with g^ 1 in (|2~4]) is 
expected in the heavy quark limit. In this limit emitted gluons will not change the spin of the 
heavy quarks, hence the spin of J/ ip is transferred to the virtual photon. Therefor the helicity of 
the two-pion state is zero. 

The function h(r, k) is unknown. We make an Ansatz for the r-dependence in the function, 
this Ansatz is motivated by the results used in the last section. We assume 



h(r,k) = a(k)(l-r) 2 (l + r) 2 , 

the function a(k) can be obtained by integrating h(r, k) over r, we obtain: 

15vr 



a(k) 



Ak° 



(7r + 7r-|a s E a (0)-E a (0)|0), 



hence for Irl < 1, 



h{r, k) 



15vr 



1 - r) 2 (l + r) 2 ■ (7r + 7r-|a s E a (0) • E a (0)|0). 



(26) 



(27) 



(2f 



The matrix element of local chromoelectric fields (7r + 7r~|a s E a (0) • E a (0)|0), which appears in the 
decay amplitude of 1 J r ' — > J/ipir + 7i~ in the QCD multipole expansion method [|T5|, [16], |l^, |18| , can 
be written in our notations as [T7J 



(7r + 7r-|a s E a (0) -E a (0)|0) 



2tt 

Y 

2tt 
~9 



(7r + 7r-\em--a s ^)M^^)(k 



2ml 



mt 



P (cos 9 V ) 



+^a s (/i)M 2 G (/i)|k| 2 /3 2 P 2 (cos 9 V 



(29) 



where 9^ v is the total energy-momentum tensor of QCD, 9® is the gluonic component of it, M 2 G (/t) 
is determined by gluon fragmentation into one pion as before. In |L7], including (9(m 2 ) corrections, 
(7r + 7r _ |^|0) = q 2 + 2m 2 is obtained from some general considerations. This coincides with the 
result of chiral perturbation theory at leading order of chiral expansion. Since the kinematic region 
we considered is only part of the whole phase space and m n7T is not very near 7r + 7r~ threshold, 
we expect that the correction from next-to-leading order of chiral perturbation theory to be 
important, so we use the expression derived from chiral perturbation theory at next-to-leading 



order for (tt+tt^ 1 1 0> , i.e. p 



7T 7T 



|0£|O> = {ml, + 2m 2 )/ (m_) + b e m 



4 

7T7T 



(30) 
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Figure 8: The differential decay branching ratio of J/ip, d BR(J/ip, S)/dm 7T7T as a function of 
m nn in unit of 10~ 5 GeV _1 with the cuts. The dashed line denotes the distribution by using the 
leading order result of chiral perturbation theory for (7r + 7r~|6^|0), while the solid line denotes the 
distribution by adding one- loop correction to the matrix element. 



with b e = 2.7GeV~ 2 . 

With these results, we are able to predict the shape of the differential decay branching ratio 
numerically. We use the cuts < |k| < j^M^ and 2m w < m nn < 0.7GeV to make the pions 
to be soft. In Fig. || The differential decay branching ratio of J/ip, d BR(J/ip, S)/dm. Kn as a 
function of m nn in unit of 10 _5 GeV _1 is shown. We use = 0.7 and M^(^) = 0.5 as used 

in 1 17]. The solid line denotes the distribution by using (30), while the dashed line denotes the 
distribution by using (7r + 7r~|6^|0) at the leading order of chiral perturbation theory. Integrating 
over 2m n < m n7T < 0.70GeV, the decay branching ratios for J/tp in the considered kinematic 
region are 1.8 x 10~ 5 and 3.5 x 10 -5 , by using the result at leading- and next-to-leading order 
of chiral perturbation theory respectively, indicating the importance of the next-to-leading order 
chiral corrections to the matrix element of the QCD total energy-momentum tensor. It should be 
noted all our numerical results are insensitive to the values of Mtfifi) and by varying the 

value of M^(/i) from to its asymptotic value (|22]) , all numerical results are changed less than 
20%. Since appears always with M| ? (/i) in the form a s (fi) ■ M^(/z), the same is also true for 

Our results indicate that this decay mode and the m nn distribution are observable at BES 
II and the proposed BES III and CLEO-C. Experiment study of the decay can test our model for 
h(r, k) or extract it. This will provide information how gluons are converted into two soft pions. 

Although we have made the Ansatz for the function h(r,k) in (26), where the shape as a 
function of r is fixed, and the parameter a(k) is just a normalization factor determined by the 
matrix element (7r + 7r~|a s E a (0) ■ E a (0)|0), but we can expect that our results for the branching 
ratio will be not changed dramatically with a change of the shape, because the normalization 
factor is fixed. 

The corresponding m nn distribution for T(IS'), referring as d BR ( T, S)/dm 7T . K , in unit of 
lO^GeV -1 < |k| < ^M T and 2m^ < m W7T < 0.7GeV is shown in Fig. |. The solid line 



12 



O 

0.3 



0.4 0.5 

m ^ ( GeV ) 



0.6 



0.7 



Figure 9: The differential decay branching ratio of T(1S), referring as d BR(T, S)/dm 7T7T , as a 
function of m nn in unit of 10~ 6 GeV _1 . The dashed line denotes the distribution by using the 
leading order result of chiral perturbation theory for (7r + 7r~|6^|0), while the solid line denotes the 
distribution by adding one- loop correction to the matrix element. 



denotes the distribution using next-to-leading order chiral perturbative theory to determine the 
matrix element (7r + 7r~|#^|0), while the dashed line denotes the distribution using leading order 
chiral perturbative theory for this matrix element. Integrating over 2m w < m W7r < 0.70GeV, the 
decay branching ratios for T(1S) in the considered kinematic region are 1.5 x 10~ 6 and 3.5 x 10~ 6 ,by 
using the result at leading- and next-to-leading order of chiral perturbation theory respectively. 
With the numerical results the decay mode may be difficult to be observed even at CLEO-C. 
However, we can learn from comparing Figs. |] and |] that when the phase space becomes larger, 
the next-to-leading order chiral corrections to the matrix element (7r + 7r~|6^|0) change the shape 
of m nn distribution dramatically. 



4 Decays of B c into a lepton pair combined with two pions 

The observation of the meson B c via the decay mode Bf — > J ftp + i^v has been reported recently 
by the Collider Detector at Fermilab (CDF) Collaboration (T^j. The B£ meson is the lowest-mass 
bound state containing a charm quark and a bottom antiquark. It has nonzero flavor and can 
decay only via weak interaction. Hence it has a very long lifetime, t(B+) = 0.46t^(stat.) ± 
0.03(syst.)ps. It will offer a new window for study the weak decay mechanism of heavy flavors 
and test various nonperturbative models for bound states. The leptonic decay of B c to one heavy 
meson has been studied in various models^], |2T|. In this section we study the leptonic decay of 
B+ into two pions. The first part of this section is devoted to leptonic decay of B^ into two hard 
pions, the decay into two soft pions is studied in the second part. 
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4.1 The leptonic decays of combined with two hard pions 

We study this exclusive decay in the rest frame of B+: 

B+(P) - l + ( Pl ) + u t (p 2 ) + tt + (K + ) + ir-(K-), (31) 

where I — e, /i, the momenta are indicated in the brackets. We study the decay in the region where 
the two-pion state has a small invariant mass and has a large total momentum. Similarly as in 
Sect. 2 the decay amplitude can be factorized, in which the nonperturbative effect related to B+ 
meson is represented by a NRQCD matrix element, and that related to the two pions is represented 
by the same distribution amplitude of two gluons in the isoscalar pion pair <3> G (a;i, (, m nn ) which 
is defined in (13). The S-matrix element for the decay is 

(f\S\i) = i^V bc L, ■ J d*ze^(n+n-\b(z)^(l - ^)c(z)\B+), (32) 

where V bc is the Cabibbo-Kobayashi-Maskawa matrix element, c(z) and b(z) is the Dirac field for 
c— quark and for 5-quark respectively, q = p\ + p 2 and 

L m = u(P2)7m(1-7 5 Mpi), (33) 

u(p 2 ) and v(pi) are the spinors of the leptons. Using the method in Sec. 2, keeping leading terms 
in heavy quark expansion and in velocity expansion, we have 

(f\S\i) = ^V bc gl{2n)U\P - k - q)L» ■ <0|x^ c |5+> 



, 1 j_ ® G (xu(,m™) 

X 



Jo 



Xi(l - Xi) 



2^t 1 ' H^ lUl (P, k, k±) 



(34) 



where xl(i>c) is the NRQCD field for b(c) quark, H^ lUl (P, k, k±) is the hard part of the decay 
amplitude and can be calculated perturbatively. We obtain: 

-dp* ■ H^ 1U1 (P, k, k ± ) = , (35) 

where = for mi = is used. The differential decay width can be written as 

1 7T 2 \k\ 

dr = . _ . JU . # El M l W<K^<"sMcos0,#, (36) 

where P' t — 1 — mf/q 2 is the velocity of l + in the center mass frame of Z + z//. 

To present numerical predictions, the NRQCD matrix element {Q\x\i> c \B^) should be known. 
It is related to B c decay constant fs c via 

mx\^ c \Bt)\ 2 = \f Bc M Ba (37) 
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with fsc = 480MeV taken from ||20|| . Other parameters take the following values: M^ c = 
6.4GeV, \V bc \ = 4.0 x W~ 2 ,G F = 1.166 x 10- 5 GeV" 2 , a s (M Bc ) = 0.24. We use the cuts: k + > 
10 k~, k° + k> 2.0GeV and 2m w < m wn < 0.70GeV. 

With these parameters and cuts we can predict the differential decay branching ratio in the 
considered region. The m nir distribution of semileptonic decay, d BR(5 C , H)/dm 7Tn in unit of 
10 _7 GeV _1 with the cuts is shown in Fig. [K| the q 2 distribution is presented in Fig. |TT[ Since the 
absolute squared matrix element of decay in this region is almost independent of q 2 , the shape 
of the q 2 distribution is determined mainly by the phase space factors. The decay branching ratio 
is 5.1 x 10~ 7 ; belonging to it is the s-wave contribution is 4.7 x 10 -7 . The estimated branching 
ratio shows that the decay mode in this region will be not observable even at the Large Hadron 
Collier (LHC). 



4.2 The leptonic decay of combined with two soft pions 

In this subsection, we study the leptonic decay of B+ combined with two soft pions. We use the 
same notation for momenta as before. With the method in [|TU| it is straightforward to obtain the 
S-matrix for the decay: 



(f\S\ih 



3V2 



H c (2vr) 4 5\P-k-q)(Q\ X \i, c \Bt) 



L„ ■ P» 



(38) 



m^rri^k ) 2 

where T n7T (k) is defined in (P4T), with our model for h(r, k) given in Sect. 3, it can be expressed as 



T n7T {k) 



5tt 



(7r + 7r-|a s E a (0)-E a (0)|0). 



(39) 
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Figure 11: The differential decay branching ratio of B£, d BK(B C , H) jdq 2 as a function of q 2 in 
unit of 10 -7 GeV~ 2 with the cuts. 



With this S-matrix element, it is straightforward to obtain the m n7T distribution of decay to 



two soft pions, which is shown in Fig. [l2|. The cuts are used: < |k| < jqMb c and Im^ < m n7T < 



0.70GeV. In Fig. [12], the solid line represents the m nn distribution by using next-to- leading order 
chiral perturbative theory to determine the matrix element (7r + 7r~|6^|0), while the dashed line 
denotes the distribution by using leading order chiral perturbative theory for this matrix element. 
Integrating over 2m 7r < m wn < 0.70GeV, the decay branching ratios for B+ in the kinematic 
region are 3.6 x 10~ 7 and 1.5 x 10~ 7 , by using the result at leading- and next-to-leading order of 
chiral perturbation theory respectively. The numerical results show that the decay mode is not 
observable even at LHC But we can learn from Figs. ^, [T2|, and |^ that when the phase spaces 
become larger, the next-to-leading order chiral corrections to the matrix element (7r + 7r~|#^|0) 
change the shape of m nn distribution dramatically. 



5 Summary 

In this paper we have studied the exclusive decay of J/ip, T, and B c into a lepton pair combined 
with two pions, where the two pions can be soft or hard with a small invariant mass. In both 
cases the decay amplitude can be factorized, in which the nonperturbative effect related to the 
heavy meson is represented by a NRQCD matrix element, and that related to the two pions is 
represented by a distribution amplitude of two gluons in the isoscalar pion pair in the case with 
hard pions, and by a correlator of chromoelectric fields in the case with soft pions. With suitable 
models for gluon conversion into soft or hard pions we are able to predict branching ratios and 
different distributions. 

Our numerical results show that the leptonic decay of J/ip combined with two hard pions or 
with two soft pions can be observed at BES II and at the proposed BES III and CLEO-C, while 
the other decays have a too small branching ratio to be observed. If the decays of J/ip are observed 
in experiment, it will provide information how gluons, which are fundamental degrees of freedom 
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Figure 12: The distribution of B~£ semileptonic decay to pion pair, referring as 

d B'R(J/ip,S)/dm 7T7T in unit of 10 _7 GeV _1 with the cuts. The dashed line denotes the distri- 
bution by using the leading order result of chiral perturbation theory for (7r + 7r~|#^|0), while the 
solid line denotes the distribution by adding one-loop correction to the matrix element. 

in QCD, are converted into observed pions. 
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